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Masses of the ground, orbitally and radially excited states of quark-antiquark 
mesons composed from the light (u, d, s) quarks are calculated within the framework 
of the relativistic quark model based on the quasipotential approach. The relativistic 
treatment of the light quark dynamics results in mass spectra which agree well with 
available experimental data for the masses of the most well-established states. The 
Regge trajectories for angular and radial excitations are constructed, and their lin- 
earity, parallelism and equidistance are verified. The assignment of experimentally 
observed light mesons to particular Regge trajectories is based on their masses and 
quantum numbers. 
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I. INTRODUCTION 

Last years an extensive analysis of the data on highly excited light non-strange meson 
states up to a mass of 2400 MeV collected by the Crystal Barrel experiment at LEAR 
(CERN) has been published [H, 0]. The classification of these new data requires a better 
theoretical understanding of light meson mass spectra. The aim of this paper is to apply 
the relativistic quark model which proved to be successful in studying various properties of 
heavy hadrons to the calculation of the masses of the radially and orbitally excited light 
meson states. All main assumptions and fixed values of model parameters are preserved 
in the present investigation. Light quarks are treated fully relativistically without the v/c 
expansion. Various non-strange and strange meson states with masses up to 2500 MeV 
are considered. This is especially important, since light exotic states (such as tetraquarks, 
glueballs, hybrids) predicted by quantum chromodynamics (QCD) are expected to have 
masses in this range [1, 0, H, |6| . The experimental data show that a large degeneracy emerges 
in the spectra of the orbitally and radially excited resonances. It is argued [7| that the states 
of the same spin with different isospins and opposite parities are approximately degenerate in 
the interval 1700-2400 MeV. An intensive debate is going on now in the literature whether the 
chiral symmetry is restored for highly excited states (see e.g. 0, @| and references therein). 
Various phenomenological and theoretical arguments, such as quasiclassical considerations, 
AdS/QCD etc. are used. 

A vast literature on the light meson spectroscopy is available. Different attempts to 
study light mesons on the basis of the relativized quark model the Dyson- Schwinger and 
Bethe-Salpeter equations 10, 11, 12], the Tamm-Dancoff method 3], chiral quark models 
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with spontaneous s ym metry breaking like the Nambu-Jona-Lasinio model IJ], finite-energy 



sum rules in QCD [l5|, lattice QCD [16|, AdS/QCD models (l7|, [l8|, etc. were undertaken. 
Therefore we mostly refer to the recent reviews where the references to earlier review and 
original papers can be found. 



In Refs. \iM, 12Q[ we studied the masses of ground and radially excited states of light 
mesons on the basis of the three-dimensional relativistic wave equation with the QCD- 
motivated potential. In this analysis we took into account the highly relativistic dynamics 
of light quarks and carried out all calculations without either the v/c oi l/m^ expansions. 
We also used the expression for the QCD coupling constant which exhibits freezing at 
small values of the momentum transfer. Good overall agreement of the obtained predictions 
and experimental data was found. The consistent relativistic treatment of the light quark 
dynamics resulted in a nonlinear dependence of the bound state equation on the meson 
mass which allowed to get the correct values of the pion and kaon masses in the model 
with the explicitly broken chiral symmetry. The obtained wave functions of the pion and 
kaon were successfully applied for the relativistic calculation of their decay constants and 
electromagnetic form factors [20|. Recently, in the framework of the same approach we 
calculated masses of the ground-state light tetraquarks using the diquark-antidiquark picture 
(2l| . It was found that scalar mesons with masses below 1 GeV agree well with the light- 
tetraquark interpretation. Indeed, it explains naturally the peculiar inverted pattern of the 
mass ordering of the lightest scalar fiavour SU{3) nonet. Here we investigate the Regge 
trajectories both in (M^,J) and (M^, n^) planes (M is the mass, J is the spin and rir is the 
radial quantum number of the meson state), check their linearity and equidistance which 
follow from experimental data [l^]. The assignment of experimentally observed mesons to 
particular Regge trajectories is proposed. 

The paper is organized as follows. In Section [TTl we describe the relativistic quark model, 
giving its main assumptions and parameters which were fixed in previous considerations. 
The relativistic quasipotential of the light quark-antiquark interaction in the meson is con- 
structed in Sec IIIII The procedure which makes this potential local and avoids the arising 
fictitious singularities is described in detail. In Sec. II VI the obtained results for the masses of 
orbital and radial excitations of light mesons are presented and compared with available ex- 
perimental data. Finally, the Regge trajectories are constructed. Their linearity, parallelism 
and equidistance is verified, and the slopes of different trajectories are compared. Section IVj 
contains our conclusions. 



II. RELATIVISTIC QUARK MODEL 

In the relativistic quark model based on the quasipotential approach a meson is described 
by the wave function of the bound quark-antiquark state, which satisfies the quasipotential 



equation of the Schrodinger type [23 



\ 2iiR 2iir) J (27r) 

where the relativistic reduced mass is 



^Af(p) = / 7^^(P' )M/M(q), (1) 



_ ^1^2 _ - {ml - mlf 
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and El, E2 are given by 

— ml + ml — ml + ml . , 

= 2M ' = 2M ■ 

Here M = Ei + E2 is the meson mass, mi^2 are the quark masses, and p is their relative 
momentum. In the center-of-mass system the relative momentum squared on mass shell 
reads 

u2(,., _ [M' - (mi + m2)^][M^ - (mi - m^f] , . 

The kernel ^(p, q; M) in Eq. ([1]) is the quasipotential operator of the quark- ant iquark in- 
teraction. It is constructed with the help of the off-mass-shell scattering amplitude, projected 
onto the positive energy states. Constructing the quasipotential of the quark-antiquark in- 
teraction, we have assumed that the effective interaction is the sum of the usual one-gluon 
exchange term with the mixture of long-range vector and scalar linear confining potentials, 
where the vector confining potential contains the Pauli interaction. The quasipotential is 
then defined by 

V{v. q; M) = Mi(p)m2(-p)V(p, q; M)Mi(g)M2(-g), (5) 

with 

V(p, q; M) = ^a.D^.(k)7r72^ + V^^M^^i"^ 2,, + Ko„f (k), 
where is the QCD coupling constant, D^^ is the gluon propagator in the Coulomb gauge 

D°°(k) = -g, D«^(k) = -ij (^5^^- - ^) , D°' = D'° = 0, (6) 

and k = p — q; 7^ and u{p) are the Dirac matrices and spinors 

u\v) ^ 



e(p) + m 1 


o-p 1 


Hv) \ 


V e{p) + m / 



X\ (7) 



with e{p) = ^pp^^^Prr? . The effective long-range vector vertex is given by 

^,(k)=7M + £^M.fc^ (8) 

where k is the Pauli interaction constant characterizing the anomalous chromomagnetic 
moment of quarks. Vector and scalar confining potentials in the nonrelativistic limit reduce 
to 

V;rnf(0 = (l-e)(Ar + 5), 

Fi„i(r) = e(Ar + 5), (9) 

reproducing 

Konf (r) = Konf (0 + Konf (0 = + 5, (lO) 

where e is the mixing coefficient. 

All the model parameters have the same values as in our previous papers 23|, |2J]. The 
light constituent quark masses = = 0.33 GeV, = 0.5 GeV and the parameters 
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of the linear potential A = 0.18 GeV^ and B = —0.3 GeV have the usual values of quark 
models. The value of the mixing coefficient of vector and scalar confining potentials e = —1 
has been determined from the consideration of charmonium radiative decays [2^. Finally, 
the universal Pauli interaction constant k = —1 has been fixed from the analysis of the fine 
splitting of heavy quarkonia ^Pj- states j23| . In this case, the long-range chromomagnetic 
interaction of quarks, which is proportional to (1 + k), vanishes in accordance with the 
fiux-tube model. 



III. QUASIPOTENTIAL OF THE LIGHT QUARK- ANTIQU ARK 

INTERACTION 

The quasipotential (jS]) can be used for arbitrary quark masses. The substitution of the 
Dirac spinors into ([5]) results in an extremely nonlocal potential in the configuration 
space. Clearly, it is very hard to deal with such potentials without any additional transfor- 
mations. In order to simplify the relativistic qq potential, we make the following replacement 
in the Dirac spinors: 



ei,2(p) = V<2 + P'^^i,2 (11) 

(see the discussion of this point in [l^, [23]). This substitution makes the Fourier transforma- 
tion of the potential ([5]) local. Calculating the potential, we keep only operators quadratic in 
the momentum acting on Vcoui, Korrf replace — >• i?^2 ~ '"^i 2 higher order operators 
in accord with Eq. (ITT!) preserving the symmetry under the (1 2) exchange. It is necessary 
to point out that such substitutions lead to the quark-antiquark potential which commutes 
with operators of the total angular momentum and the orbital angular momentum. There- 
fore J and L are good quantum numbers, as in the nonrelativistic approach. However the 
nonlinear dependence of the interaction potential on the meson mass effectively takes into 
account the relativistic character of the light quark interaction. Note that the global features 
of highly excited light mesons can be well understood in terms of the relativistic relations 



involving J as well as nonrelativistic relations involving L |17l . 1251 . 

The substitution ffTTj) works well for the confining part of the potential. However, it leads 
to fictitious singularities and 6^{r) at the origin arising from the one-gluon exchange 
part ( A Vcoui ('")), which is absent in the initial potential. Note that these singularities are 
not important if they are treated perturbatively. Since we are not using the expansion in 
v/c and are solving the quasipotential equation with the complete relativistic potential, an 
additional analysis is required. Such singular contributions emerge, e.g., from the following 
terms 



[ei(g)ei(p)e.(g)e,(p)]V2^^-'(k^)' (12) 

if we simply replace ei,2 -£'1.2- However, the Fourier transforms of expressions (JT2l) are 
less singular at r — 0. To avoid such fictitious singularities we note that if the binding 
effects are taken into account, it is necessary to replace ei_2 -Ei,2 — Vi,2V, where V is 
the quark interaction potential and 771,2 = ""^2,1/(^1 + ^2)- At small distances r —>■ 0, the 
Coulomb singularity in V dominates and gives the correct asymptotic behaviour. Therefore, 
we replace ei,2 —>■ Ei,2 — ''?i,2Vcoui in the Fourier transforms of terms (fT2l) (cf. [2^). We 
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used a similar regularization of singularities in the analysis of heavy-light meson spectra 
(23 |. Finally, we ignore the annihilation terms in the quark potential since they contribute 
only in the isoscalar channels and are suppressed in the ss vector channel. 
The resulting qq potential then reads 



V{r) = Vsi{r) + Vsj,{r), 
where the spin-independent potential has the form 

{Ef -ml + El - 



(13) 



Vsiir) 



Vcoul(r) + K;onf(r) + 
1 



+ 



77111712 

El + mi 
El 



1 + il + K) 



4(Ei + mi)(E2 + m2) [E1E2 
(El + mi)(E2 + m2) 



Vcoui(r) 



(I + k) 



E1E2 



+ 



El + 7712 
Eo 



mim2 



-K 



conf \ 



1 

+ - 



4 \Ei{Ei + mi) 
1 

'4 



mi(_E'i+mi) 77i2{E2 + 7712) 
\El -ml + El - mi) 



E2{E2 + m2) 
- (1 + /^ 



(2) 



CouU 



/ 1 



8mi77i2{Ei + mi)(i?2 + "^2) 
and the spin-dependent potential is given by 



\Ei7ni 
1 L2 



E27n2 



AKrnf(r) 



E1E2 2r 



^Coul(^), 



(14) 



^sd('^) = cti LSi + a2 LS2 + b 



-SiS2 + ^(Sir)(S2r) 



+ cSiS2 + ci(LSi)(LS2), (15) 



ai 



2E1E2 



X 



/Ei 



(2. 

- mi 



2mi 



i^i + mi / r _fci + mi r 

^Ei + mi\ 2 ( ^^fEi - 



2m2 



AE1E2 {El 
fE2 ^, 



02 
h 



a + 

X ( — - 2(1 + 
Vm2 

ai(1^2), 

^ ^^Coui(^: 



1 + ^^) 

Ei + mi\ 2 y ( El - mi E2-m2\ 1 / 

2mi J r V-c/2 + ^2 -C/i+mi/r 

-^TTT^ T V-y'coAr) + -Konf W + f — - 2(1 + «)^^^ 

mi)(£;2 + m2) [r ' r ^ Vmi ^ ^ 2mi 

^2 + m2\ / 



2m2 



mi' 



3E1E2 ^ 
f E2 — 

X 



m2 



2m2 



-(1 + 



K 



3E1E2 



X 



AVconiir) + 



-(1 + 



^coui(r) + 

E2 + 7712" 
2m2 
^i?i — mi 



f El — mi 



2m 

I „ ^on 



-(l + «:)^ 
^ ^ 9m 



mi ' 



2mi 



/'E2- 77l2 



K 



\ 2m2 
^1^2(^1 +mi)(E2 + m2) 



(1 + k) 

2mi 2mi 

Ei + m^ 



1 



2m2 

1 



^coui(0 - v;^'oui(0 + \vu{r) - KLi{r) 



(16) 
(17) 

(18) 

(19) 

(20) 
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with 



_4a^ 
3 r ' 

Vboui(^) 



l + Vi 



4a, 1 
3Eir 



l + Vi 



a. 



3 Ei + nii r 



1 + Vi 



3Eir 



1 + V2 



3E2r 
V{r 



Vl,2 



1,2), 



m2A 



mi + m2 



3Eir 



4 



3 El + mi r 



1+V2 



4a, 1 
3E2r 



1+V2 



1 



3E2 + m2r 



(21) 



Here we put a, = a. 



with yUi2 = 2mim2/(mi + m2). We adopt for a,(/i^) the simplest 



model with freezing |27|, namely 

a,(/i^) = - 



An 



A2 



2 ' 



11 



(22) 



where the background mass is Mb = 2.2AV A = 0.95 GeV |27|, and A = 413 MeV was fixed 
from fitting the p mass. ^ We put the number of flavours nj = 2 for vr, p, K, K* and nj = 3 
for 0. As a result we obtain as(/i^^) = 0.730, as(/i^g) = 0.711 and as{p1s) = 0.731. Note 
that the other popular parametrisation of a, with freezing 28| gives close values. 



IV. RESULTS AND DISCUSSION 

The calculated masses of light unfiavoured and strange mesons are given in Tables HI and HTl 
They are confronted with available experimental data from PDG Particle Listings including 
data from the "Further States" Section [2^. We find good agreement of our predictions 
with data. Most of the well-established state masses are reproduced in our model. 

We do not consider the mixing of states in the isoscalar sector. Therefore the predictions 
in Table [T] are given for the pure qq and ss states. Such mixing is mostly important in 
the pseudoscalar sector. We follow the rj — rj' mixing scheme proposed in Ref. and 
take the phenomenological values of the mixing angle (p = 38° and the decay constant ratio 
y = fq/ fs = 0.81. Using our values for the mass of M^^- and the pion mass we get = 
573 MeV and M^/ = 989 MeV close to the measured masses M^^p = 547.853 ± 0.0024 MeV 
and M^f^ = 957.66 ± 0.24 MeV jl^. The experiment shows that the vector and excited 
isoscalar light mesons are almost ideally mixed and therefore can be roughly considered as 
pure qq and ss states. Indeed we find reasonable agreement of our prediction with experiment 
in the isoscalar sector. 



The definition ([^ of as can be easily matched with the as used for heavy quarkonia [23| at the scale 
about rric. 



TABLE I: Masses of excited light {q = u,d) unflavored mesons (in MeV). 







Theory 

qq 




Experiment 




Theory 


Experiment 


/ = 1 


mass 


/ = 


mass 


ss 


/ = 


mass 




0-+ 


154 


TT 


139.57 






743 








1— 


776 


P 


775.49(34) 


u 


782.65(12) 


1038 




1019.455(20) 




0++ 


1176 


do 


1474(19) 


fo 


1200-1500 


1420 


fo 


1505(6) 


l^Pi 


1++ 


1254 


ai 


1230(40) 


fl 


1281.8(6) 


1464 


fl 


1426.4(9) 


1^P2 


2++ 


1317 


0,2 


1318.3(6) 


f2 


1275.1(12) 


1529 


fi 


1525(5) 


l^Pl 


1+- 


1258 


hi 


1229.5(32) 


hi 


1170(20) 


1485 


hi 


1386(19) 




0-+ 


1292 


TT 


1300(100) 


V 


1294(4) 


1536 


V 


1476(4) 




1— 


1486 


P 


1465(25) 


u 


1400-1450 


1698 




1680(20) 




1— 


1557 




1570(70) 


u 


1670(30) 


1845 








2~ 


1661 










1908 






13^3 


3 


1714 


P-i 


1688.8(21) 




1667(4) 


1950 




1854(7) 




2-+ 


1643 




1672.4(32) 


m 


1617(5) 


1909 


m 


1842(8) 




0++ 


1679 






fo 


1724(7) 


1969 






2^ Pi 


1++ 


1742 


a\ 


1647(22) 






2016 


fl 


1971(15) 


2^P2 


2++ 


1779 


0-2 


1732(16) 


f2 


1755(10) 


2030 


f2 


2010(70) 


2^Pi 


1+- 


1721 










2024 








0-+ 


1788 


TT 


1816(14) 


V 


1756(9) 


2085 


V 


2103(50) 


3^Si 


1— 


1921 


P 


1909(31) 




1960(25) 


2119 


v 


2175(15) 


l^F2 


2++ 


1797 






h 


1815(12) 


2143 


f2 


2156(11) 




3++ 


1910 


as 


1874(105) 






2215 


fs 


2334(25) 


13P4 


4++ 


2018 




2001(10) 


h 


2018(11) 


2286 






llPg 


3+- 


1884 










2209 




2275(25) 


2^Di 


1— 


1895 


p 


1909(31) 






2258 


LO 


2290(20) 


2^D2 


2~ 


1983 


P2 


1940(40) 


U2 


1975(20) 


2323 






2^Ds 


3— 


2066 










2338 






2^D2 


2-+ 


1960 




1974(84) 


m 


2030(20) 


2321 






S^Po 


0++ 


1993 


do 


2025(30) 


fo 


1992(16) 


2364 


fo 


2314(25) 


33Pi 


1++ 


2039 


ai 


2096(123) 






2403 






3^P2 


2++ 


2048 


0'2 


2050(42) 


f2 


2001(10) 


2412 


f2 


2339(60) 


3^ Pi 




2007 


bi 


1960(35) 


hi 


1965(45) 


2398 








0-+ 


2073 


IT 


2070(35) 


V 


2010(50) 


2439 








1— 


2195 


P 


2265(40) 


u 


2205(30) 


2472 








3— 


2002 


P3 


1982(14) 


UJ3 


1945(20) 


2403 








4~ 


2122 


Pi 


2230(25) 


U4 


2250(30) 


2481 








5~ 


2264 


P5 


2300(45) 


UJ5 


2250(70) 


2559 








4-+ 


2092 










2469 






3^Di 


1— 


2168 


P 


2149(17) 






2607 






S^D2 


2— 


2241 


P2 


2225(35) 


UJ2 


2195(30) 


2667 








3— 


2309 


P3 


2300(60) 




2278(28) 


2727 






3^D2 


2-+ 


2216 


7r2 


2245(60) 


m 


2248(20) 


2662 
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TABLE I: (continued) 



n 







Theory 




Experiment 




Theory 






QQ 


/ = 1 


mass 


1 = 


mass 


ss 


2^F2 


2++ 


2091 


a2 


2100(20) 




2141(12) 


2514 


23^3 


3++ 


2191 


«3 


2070(20) 






2585 


23^4 


4++ 


2284 






h 


2320(60) 


2657 


2IF3 


3+- 


2164 


bs 


2245(50) 






2577 


43Po 


0++ 


2250 






fo 


2189(13) 


2699 


43 Pi 


1++ 


2286 


ai 


2270(50) 


fl 


2310(60) 


2729 


43P2 


2++ 


2297 


02 


2280(30) 


f2 


2297(28) 


2734 


4^ Pi 


1+- 


2264 


bi 


2240(35) 


hi 


2215(40) 


2717 


23^3 


3— 


2267 


P3 


2260(20) 


UJ3 


2255(15) 


2743 


23^4 


4— 


2375 










2819 




5— 


2472 










2894 


2IG4 


4-+ 


2344 




2250(15) 




2328(30) 


2806 




0-+ 


2385 


vr 


2360(25) 


V 


2320(15) 


2749 


53^1 


1— 


2491 










2782 


l3i74 


4++ 


2234 


04 


2237(5) 


fj 


2231.1(35) 


2634 




5++ 


2359 










2720 




6++ 


2475 




2450(130) 


k 


2465(50) 


2809 


l'H5 


5+- 


2328 










2706 



/ = 



mass 



The strange meson states (L^) with J = L, given in Table [TTl are the mixtures of 
spin-triplet (^Ll) and spin-singlet (}Ll) states: 

Kj = K(^Ll) cosif + K{^Ll) simp, 

K'j = -KCLL)smip + K{^LL)cosip, J = L = 1,2,3... (23) 

where is a mixing angle. Such mixing occurs due to the nondiagonal spin-orbit and tensor 
terms in Eq. f|T5l) . The masses of physical states were obtained by diagonalising the mixing 
terms. The found values of mixing angle if are the following: IP 43.8°, 2P 44.6°, 3P 44.8°, 
ID 44.2°, 2D 44.5°, IF 44.3°, IG 44.3°. These values show that physical Kj mesons are 
nearly equal mixtures of the corresponding spin-singlet Kj^Lj,) and spin-triplet K{^Ll) 
states which is in good agreement with the experimental data |29|. 

The scalar sector presents a special interest due to its complexity and the abundance of 
experimentally observed light states. We see from Table [T] that the masses of the lightest 
qq scalar mesons have values about 1200 MeV. This confirms the conclusion of our recent 
consideration [2H that light scalar mesons, /o(600) (cr), ^^0(800) (k), /o(980) and ao(980), 
with masses below 1 GeV should be described as light tetraquarks consisting from scalar 
diquark and antidiquark. Moreover the predicted masses of the scalar tetraquarks composed 
from axial- vector diquark and antidiquark |2l| have masses in the same range as the lowest qq 
scalar mesons. Therefore mixing between these states can occur, e.g. due to the instanton- 
induced mixing terms [3l|, [si. The obtained results for the masses indicate that ao(1450) 



should be predominantly a tetraquark state which predicted [21| mass 1480 MeV is within 
experimental error bars MaQ(i45o) = 1474 ± 19 MeV. The exotic scalar state X(1420) from 
the "Further States" Section could be its isotensor partner. On the other hand sqil^Po) 
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TABLE II: Masses of excited strange mesons (in MeV). 
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lliClOO 
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0" 
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0^ 
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1~ 
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K* 
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1^ 
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1^ 
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0+ 


1424 
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— 

3 
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1412 
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2 
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Ko 
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1^ 
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2 
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0~ 


1538 






3^Po 


0^ 
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1" 
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K* 
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2+ 
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1^ 
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K* 
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3Pi 
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2200 
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6 
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1789 


^3 
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3Pi 
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2 


1824 


Ko 


1816(13) 


1 


3 


2207 






11^2 


2 


1 lUy 


Kn 


1 773l'8'> 


1 (-^5 


r — 



2356 


Ki 


2389('94'l 


o3 7~> 


n+ 

0^ 
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interpretation is favored for Kq{1A?)Q) (see Table HTl) . This picture naturally explains the 
experimentally observed proximity of masses of the unflavoured ao(1450) and /o(1500) with 
the strange ii^^(1430). Therefore one could expect an additional isovector predominantly qq 
state ao with the mass about 1200 MeV, though it was not observed in several experiments. 

It was noted long ago that the light meson Regge trajectories are almost linear in (J, M^) 
and (?T,r,M^) planes. 

a) The (J, M^) Regge trajectory: 

J = aM^ + ao; (24) 

b) The {rir, M^) Regge trajectory: 

nr = n-l = (3M'^ + (3o, (25) 

where a, /3 are the slopes and ao, /?o are intercepts. The relations (p^ and (!25|) arise in 
most models of quark confinement, but with different values of the slopes. For example, the 
QCD string with two light quarks at the ends gives the slopes (33| : 
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5 : 
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J 

FIG. 1: Parent and daughter {J,M'^) Regge trajectories for isovector light mesons with natural 
parity (p). Diamonds are predicted masses. Available experimental data are given by dots with 
error bars and particle names. is in GeV^. 




where a is the string tension which is equal to the slope of the linear confining potential A 
in Eq. (fTO|). 

On the other hand, the quasiclassical picture for a light meson, described by the massless 
Salpeter equation with a linear confining potential: 



(2p + Ar)'ip = Mip, 



(27) 



gives for the Regge slopes 



a 



1 

8A' 



in A' 



(28) 



implying that 



a/p = n/2. 



(29) 

In Figs. [THDand Figs. [MS] we plot the Regge trajectories in the (J, M^) plane for mesons 
with natural (P = (—1)"^) and unnatural (P = (—1)'^^^) parity, respectively. The Regge 
trajectories in the (n^, M^) plane are presented in Figs. fT0lfT2l The masses calculated in 
our model are shown by diamonds. Available experimental data are given by dots with 
error bars and corresponding meson names. Straight lines were obtained by a fit of the 
calculated values. The fitted slopes and intercepts of the Regge trajectories are given in 
Tables HII] and [IVl We see that the calculated light meson masses fit nicely to the linear 
trajectories in both planes. These trajectories are almost parallel and equidistant. 

It is important to note that the quality of fitting the vr meson Regge trajectories both in 
(J,M2) and {rir.M^) planes is significantly improved if the ground state vr is excluded from 
the fit (the is reduced by more than an order of magnitude and becomes compatible with 
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TABLE III: Fitted parameters of the {J,M'^) parent and daughter Regge trajectories for light 



mesons with natural and unnatural parity {q = u,d). 



Trajectory 


natural parity 


unnatural parity 




Ct 1 vjrc V j 


OLq 


IX 1 VJTC V 1 




QQ 


H 




TT 




parent 


fi QQ7 _|_ fi nno 

W.OO 1 —1— \J.\J\J\J 


n AKR + n 018 




W.W.^(J —I— W.UtJT: 


daughter 1 


1.009 ± 0.019 


-1.232 lb 0.074 


1.031 zb 0.063 


-1.846 zb 0.217 


daughter 2 


1.144 lb 0.113 


—3 092 zb 540 


1 171 zb 009 


—3 737 zb 042 


qq 










parent 


L.LZiO ZIZ U.UOO 


_i acvi -U n 1 rizi 


j-.ui^ zn u.uou 


— U.UOO zc U.l^U 


daughter 1 


1.291 ± 0.003 


-3.640 ± 0.011 


1.148 zb 0.012 


-2.497 ± 0.050 


daughter 2 


1.336 ± 0.022 


-5.300 ± 0.102 


1.154 zb 0.014 


-3.798 ± 0.007 


qs 


K* 




K 




parent 


0.839 ± 0.004 


0.318 ±0.012 


0.780 zb 0.022t 


-0.197 zb0.036t 


daughter 


0.942 lb 0.046 


-1.532 zb 0.209 


0.964 zb 0.072 


-2.240 zb 0.296 


ss 










parent 


0.728 zb 0.011 


0.234 zb 0.034 


0.715 zb 0.023 


-0.444 zb 0.068 


daughter 1 


0.721 ± 0.089 


-1.072 zb 0.047 


0.718 ± 0.032 


-1.786 ±0.157 


daughter 2 


0.684 ± 0.039 


-2.047 zb 0.226 


0.729 ± 0.010 


-3.174 ±0.057 




* fit without tt: 


a = (1.053 ± 0.059) GeV' 


-2 , ao = -0.725 ± 0.170 




^ fit without K: 


a = (0.846 ± 0.013) GeV" 


-2 , ao = -0.431 ± 0.042 



TABLE IV: Fitted parameters of the {nr,M^) Regge trajectories for light mesons. 



Meson 


P (GeV-2) 


Po 


Meson 


p (GeV-2) 


Po 








ss 






TT 


0.679 ± 0.023* 


-0.018 ±0.014* 


Vss 


0.559 ± 0.009 


-0.315 ± 0.026 




0.700 ± 0.023 


-0.451 ±0.060 




0.597 ±0.009 


-0.662 ±0.031 


do 


0.830 ± 0.032 


-1.214 ±0.109 


fo 


0.566 ± 0.009 


-1.156 ±0.039 


ai 


0.840 ± 0.037 


-1.401 ±0.134 


fi 


0.561 ±0.013 


-1.224 ±0.058 


h 


0.863 ± 0.030 


-1.431 ±0.106 


hi 


0.575 ± 0.015 


-1.292 ±0.066 


a2(3P2) 


0.867 ± 0.036 


-1.585 ±0.134 


f2 


0.581 ± 0.007 


-1.370 ±0.031 




0.894 ± 0.013 


-2.182 ± 0.050 










0.916 ± 0.032 


-2.514 ±0.134 










0.874 ± 0.041 


-2.623 ±0.189 








02(^^2) 


0.891 ±0.010 


-2.881 ±0.043 








as 


0.890 ±0.014 


-3.254 ±0.066 








^3 


0.906 ±0.015 


-3.225 ± 0.071 








04 


0.899 ± 0.016 


-3.672 ± 0.084 









* fit without tt: P = (0.750 ± 0.032) GeV-^, Pq = -0.287 ± 0.109 
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FIG. 2: Same as in Fig. [T] for isoscalar light qq mesons with natmal parity {uo) 
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FIG. 3: Same as in Fig. [T] for isovector hght qq mesons with natural parity (ao)- 

the values for other trajectories). In the kaon case omitting the ground state also improves 
the fit but not so dramatically as for the pion. The corresponding trajectories are shown in 
Figs. O El and [10] by dashed lines, the fitted values of the slopes and intercepts are given in 
the footnotes to Tables IIIII and IIV[ This indicates the special role of the pion originating 
from the chiral symmetry breaking. 
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FIG. 4: Same as in Fig. [T] for isodublet light mesons with natural parity {K* 
M2 




J 

FIG. 5: Same as in Fig. [1] for isoscalar light ss mesons with natural parity (ip). 

It can be seen in Figs. [T]and[3]that p(1700) and ao(1450) do not lie on the corresponding 
Regge trajectories. This further confirms our previous conclusion that ao(1450) should be 
predominantly a tetraquark state and suggests the possible exotic nature of p(1700). 

From the comparison of the slopes in Tables UTTl IIVI we see that the a values are system- 
atically larger than the f3 ones. The ratio of their mean values is about 1.3 both for the light 
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M2 




FIG. 6: Same as in Fig. [T] for isovector light mesons with unnatural parity (vr). Dashed line 
corresponds to the Regge trajectory, fitted without n. 



qq isovector and ss mesons. Such ratio is lower than predictions of the QCD string model 
( !26|) and massless Salpeter equation (!29|) . The mean value of the slope (3 ~ 0.85 GeV~^ for 
isovector mesons is about two times larger than the result of the above mentioned models 
/3 = 1/{AitA) ^ 0.44 GeV"^ for A = 0.18 GeV^ used in our approach. 

The assignment of the experimentally observed states to the corresponding Regge tra- 
jectories in our model based on their masses and J^*" values (see Figs. [Tl fT2|) is slightly 
different from the previous phenomeno logical analysis j22|, 25 1 based on the equal values for 
the slopes a and /3. However the number of states, for which such correspondence is found, 
is approximately the same. Future experimental data will shed further light on this issue. 



V. CONCLUSIONS 

The mass spectra of light quark-antiquark mesons were calculated in the framework of 
the QCD-motivated relativistic quark model. All considerations were carried out without 
application of the unjustified nonrelativistic v/c expansion. Such approach leads to the 
nonlocal quasipotential of the relativistic quark-antiquark interaction. To make it local, the 
substitution (fTTj) was used. As a result the relativistic local quasipotential was obtained 
which depends on the mass of the meson in a complicated nonlinear way. Such dependence 
allowed us to get masses of the vr and K mesons in agreement with experimental data in the 
considered model, where the chiral symmetry is explicitly broken by the constituent quark 
masses. It was found that the lightest scalar (1^-Po) states have masses above 1 GeV. This 
confirms our previous conclusion jlH that /o(600) (a), Kq{800) (k), /o(980) and ao(980) 
could be the diquark-antidiquark tetraquark states. It was found that the calculated masses 
of light mesons reproduce the linear Regge trajectories both in the (J, M^) and (nr,M^) 
planes. Their slopes and intercepts were determined. The slope of the orbital excitations 
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FIG. 7: Same as in Fig. [T] for isovector light qq mesons with unnatural parity (ai). 
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FIG. 8: Same as in Fig. [T] for isodublet light mesons with unnatural parity {K\ Dashed line 
corresponds to the Regge trajectory, fitted without K. 

a was found to be in average 1.3 times larger than the slope of the trajectories of radial 
excitations (3. This value of the ratio a/ (3 is smaller than the predictions based on the 
spinless Salpeter equation fl27|) and the QCD string model. The obtained slope (3 for the 
isovector light mesons is almost two times higher than the value predicted by the above 
models. Possible experimental candidates for the states populating the Regge trajectories 
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M2 




J 

FIG. 9: Same as in Fig. [1] for isoscalar light ss mesons with unnatural parity. The ground state 
with J = is the mixture of r] and 77' with pure ss quark content {rfss)- 



are identified in Figs. [TlfT2l Predictions for the masses of the missing states are presented 
in Tables [IflTl It is clearly seen that the chiral symmetry is not restored for highly excited 
states in our model. This should be expected since the Lorentz-scalar part of the confining 
potential explicitly breaks the chiral symmetry. Our results in some cases differ from the 
previous phenomenological prescriptions (2^ . 25| . Future experimental data can help in 
discriminating between the theoretical predictions. 
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and A. Vainshtein for support and discussions. This work was supported in part by the Rus- 
sian Science Support Foundation (V.O.G.) and the Russian Foundation for Basic Research 
(RFBR), grant No.08-02-00582 (R.N.F. and V.O.G.). 
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